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ABSTRACT 

It is now clear that there is a substantial population of primordial binaries in galactic 
globular clusters and that binary interactions arc a major influence on globular cluster 
evolution. Collisional interactions involving stars in binaries may provide a significant 
channel for the formation of various stellar exotica, such as blue stragglers, X-ray bina- 
ries and millisecond pulsars. We report on an extensive series of numerical experiments 
of binary-binary scattering, analysing the cross-section for close approach during in- 
teractions for a range of hard binary parameters of interest in globular cluster cores. 
We consider the implied rate for tidal interactions for different globular clusters and 
compare our results with previous, complementary estimates of stellar collision rates 
in globular clusters. We find that the collision rate for binary-binary encounters dom- 
inates in low density clusters if the binary fraction in the cluster is larger than 0.2 for 
wide main-sequence binaries. In dense clusters binary-single interactions dominate the 
collision rate and the core binary fraction must be ^ 0.1 per decade in semi-major axis 
or too many collisions take place compared to observations. The rates are consistent if 
binaries with semi-major axes ~ 100 AU are overabundant in low density clusters or if 
breakup and ejection substantially lowers the binary fraction in denser clusters. Given 
reasonable assumptions about fractions of binaries in the cores of low density clusters 
such as NGC 5053, we cannot account for all the observed blue stragglers by stellar 
collisions during binary encounters, suggesting a substantial fraction may be due to 
coalescence of tight primordial binaries. 
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1 INTRODUCTION 

As the evidence for the presence of primordial binaries in 
globular clusters increases, it has become clear that the con- 
tribution of binary-single star and binary-binary scattering 
to stellar collisions and other stellar binary processes must 
be significant, at least in some clusters (see reviews by Hut 
et al. 1992; Livio 1995; also, Sigurdsson & Phinney 1995, 
Davies 1995, Davies & Benz 1995, Leonard 1989, Goodman 
& Hut 1989). Of particular importance are tidal encoun- 
ters, or stellar collisions, that occur during resonances that 
develop during hard binary-single and binary-binary scat- 
terings. These may contribute significantly to the formation 
of blue stragglers (Leonard 1989, Leonard & Fahlman 1991, 
Leonard & Linnell 1992), X-ray binaries, MSPs, CVs run- 
away stars and other exotica (Sigurdsson & Phinney 1995, 
Davies 1995). 

Binary-binary scattering may be particularly impor- 
tant in low-density clusters where there may be a large 
number of primordial binaries and products of binary in- 
teractions such as blue stragglers (Hills 1975, Nemec & Har- 



ris 1987, Nemec & Cohen 1989, Leonard 1989, Mateo et al. 
1990, Bolte 1991, Hills 1992, Leonard & Linnell 1992, Bolte 
et al. 1993, Yan & Mateo 1994). Mass segregation effects in 
globular clusters will increase the binary fraction in the core 
compared to the rest of the cluster. Thus, even if the binary 
fraction in the whole cluster is low (say ~ 5%) the fraction 
in the core may be much higher (see, for example, Leonard 
1989, Hut et al. 1992, McMillan & Hut 1994, but note also 
Sigurdsson & Phinney 1995). 

Here we report the results of 100,000 numerical exper- 
iments of hard binary-binary scatterings, for a range of bi- 
nary parameters appropriate to globular cluster interactions. 
Other studies of binary-binary encounters have been car- 
ried out (Mikkola 1983, 1984a,b, Hoffer 1983, Leonard 1989, 
Leonard & Fahlman 1991, McMillan et al. 1990, 1991, Hut et 
al. 1992, Hut 1995, Rasio et al. 1995). Mikkola considered a 
range of hard and soft binary scatterings looking at the final 
state and energy transfer, while Hoffer included mostly soft 
binary encounters. Leonard's work overlaps with ours, but 
does not present a systematic survey of cumulative cross- 
sections as reported here, and our work should be consid- 



Binary-Binary Scattering 2 



ered complementary to his. McMillan, Hut and Rasio have 
so far mostly reported studies of particular sets of encoun- 
ters or encounters in particular models of clusters rather 
than surveys of cross-sections. We present a set of cumula- 
tive cross-sections for close approach during hard encoun- 
ters for a range of mass ratios and semi-major axis ratios. 
We compute the relative event rate for the various possible 
outcomes of the encounters. We also present sample cross- 
sections for the change in semi-major axis during flybys and 
compare them to the one seen in encounters between bina- 
ries and single stars. We leave a detailed discussion of the 
subsequent evolution of the systems produced in encounters, 
such as triple-star systems, to a later paper. 

The cross-sections for close approaches calculated here 
complement previous hydrodynamical calculations of the 
outcome of stellar collisions and strong tidal interactions in 
the context of hard binary encounters (Davies et ed. 1994, 
Davies & Benz 1995, Goodman & Hernquist 1991, Sigurds- 
son & Hernquist 1992). 



2 METHOD 

The initial conditions for the scatterings were set following 
the method of Hut & Bahcall (1983, see also Sigurdsson & 
Phinney 1993). With two binaries, we have additional pa- 
rameters from the relative phase of the second binary, the 
orientation of the plane of the second binary and the sec- 
ond binary mass ratio, semi-major axis and eccentricity. We 
drew the binary parameters by Monte Carlo selection uni- 
formly over the phase variables. The relative velocity at in- 
finity of the centres-of-mass of the two binaries, Vcc was 
chosen uniformly on the interval allowed. We refer to a set 
of encounters performed at fixed a;, M; and range of 

^oo as a 

"run" . A discrete set of values for the semi-major axis, ai , 02 
and masses, Mi, M2, M 3 , M4 was used for each run. The bi- 
nary eccentricities, ei,2 were zero for all encounters reported 
here; previous calculations indicate the cross-sections of in- 
terest are not sensitive to the binaries' eccentricities. We 
discuss the limitations of this assumption later in this pa- 
per. For the runs discussed here the binary parameters used 
are shown in Table 1. 

The critical velocity, v c , is the velocity for which the 
total energy of the system in the centre-of-mass frame is 
zero, is given by 



Table 1. The properties of all the runs performed. 
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where n = (Mi+M 2 )(M 3 + M4)/M T , Mr = M1+M2 + M3 + 
M4, is the binaries reduced mass, and a\,a 2 are the semi- 
major axes of the binaries containing masses Mi i2 ,M 3i 4 re- 
spectively. Note that for these simulations Mi = Ms and 
M2 = M4. As a convention we order ai > ai, Mi, 3 > Mm, 
and choose G = 1. The sampling in velocity was uniform in 
Vac/vc over the range indicated for each set of runs shown. 

We refer to encounters where Voa/v c < 1 as "hard", fol- 
lowing the nomenclature established for binary-single scat- 
terings. Hard encounters are dominated by gravitational fo- 
cusing. Treating the binaries as point masses at each binary 
centre-of-mass, for an impact parameter 6, the pericentre, 
p is given by 
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p = 



GM T 



1/2 



(2) 



For Voo/vc <l,p» b 2 v 2 a /2GMT- The impact parameter for 
each scattering is uniform in b to some maximum impact 
parameter b max . By extension of Hut & Bahcall's choice 
(1983) we set b max — Cai/voo +Dai, where C — 5, D — 0.6 
for the set of runs reported here. For Voo/v c <§C 1, 



Pmax — 



C 2 aj 
2GM T 



1 + 



C 



(3) 



Note for 01 3> 02 the maximum pericentre approach is large 
compared to 02; this is necessary as the wider binary may 
be sensitive to perturbations from the tighter binary at sev- 
eral 01, while the tighter binary will likely be only very 
weakly perturbed. Thus for ai /a 2 ^> 1 we have to sam- 
ple the scatterings to large impact parameter. For a\ = 4d2 
and 01 = 16a2 we carried out a separate run with C = 4 to 
check the cross-section for very close approaches had con- 
verged and we were sampling the strong interactions ade- 
quately. The runs with C = 5 proved adequate and results 
from those are reported here to provide a homogenous sam- 
ple. The results from the smaller impact parameter runs will 
be discussed in a later paper. 

While the simulations are scale free, the choice of masses 
and Voo/vc were made bearing in mind the physically inter- 
esting range of velocities in globular clusters, Voo ~ 10km s _1 
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and Mi — 0.5 — 1.5MQ. Of particular interest are wide 
(o» ~ 1 — 100 AU) binaries containing main-sequence stars 
near the turnoff (M, « 0.7M©), neutron stars (Mi « 1.4M©) 
and white dwarfs (Mi ~ 0.5 - 1.2M ). 



2.1 Integration scheme 

Two integration schemes were used in the calculations: a 
fourth order Runge-Kutta integration scheme with adap- 
tive step size and quality control (see Hut & Bahcall 1983, 
Sigurdsson & Phinney 1993), and a Bulirsch-Stoer variable 
step integrator with KS-chain regularisation (Aarseth 1984, 
Mikkola 1983, 1984a, b). The Runge-Kutta scheme is sim- 
ple to implement and provided a direct comparison with 
previous binary-single scatterings. However, for ai/a2 ^ 2 
the step size necessary to prevent secular drift in the total 
energy of the system becomes prohibitively small, and the 
Bulirsch-Stoer regularised scheme is an order of magnitude 
faster in integration despite the higher cost per integration 
step. The Bulirsch-Stoer regularised scheme is more com- 
plicated to implement and we relied heavily on subroutines 
provided by Sverre Aarseth. Typically the Runge-Kutta in- 
tegration for ai = 0,2 required <; 10 steps to resolve a hard 
encounter, while the Bulirsch-Stoer regularised scheme typ- 
ically required <; 10 4 steps. 

As the Bulirsch-Stoer scheme uses large integration 
steps, there is concern that it may not accurately track the 
true close approach pair separations. We note that the cross- 
sections from the set of runs carried out with the Runge- 
Kutta integration scheme agreed to within statistical error 
with the cross-sections calculated by the Bulirsch-Stoer reg- 
ularised integrator for the same sets of parameters. We also 
varied the parameter for integration tolerance by two orders 
of magnitude for one set of runs, forcing a smaller integration 
step size, and checked that the cross-sections did not change 
with the integration step size. Here we report a homogenous 
set of runs done only with the Bulirsch-Stoer regularised 
scheme. 

The code was run on a DEC 3000/400 alpha, and the 
total set of runs required about 3 weeks of cpu time. During 
each encounter the true pairwise separation between each 
pair of particles was monitored and the minimum value of 
each was stored. In addition we stored the position and ve- 
locity of all four particles at the moment of the single closest 
pair approach. This was for future analysis of the properties 
of the remnant system for particles deemed to have under- 
gone a strong tidal encounter of collision. As the simulations 
are scale free, we chose not to pick a scale and allow tidal in- 
teractions to occur during the encounter, rather we analyse 
the outcome using a "sticky particle approximation" with 
a variable scale picked for each run after the runs are com- 
pleted (Sigurdsson & Phinney 1993, Davies et aJ. 1994). The 
analysis using the sticky particle approximation will be dis- 
cussed in a later paper (in preparation). The position and 
velocities saved at closest approach also permit generation of 
initial conditions for SPH simulations by time reversing the 
integration from the point of closest approach to a suitable 
initial separation (see Sigurdsson & Hernquist 1992). 

The integration was stopped when either a maximum 
number of integration steps had been taken, 2,000,000 for 
Runge-Kutta integration, 500,000 for Bulirsch-Stoer in- 
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Figure 1. An example of an exchange during an encounter be- 
tween unequal mass binaries with different semi-major axis. 



tegration; or when two of the pair separations, nj, ex- 
ceeded some critical value R m = 1.2 x R in where R in — 
30 x max{ai,a2} was the initial binary separation. 

2.2 Possible Outcomes 

There are several possible outcomes of encounters between 
two binaries. In a flyby the two binaries may remain in- 
tact, with the components remaining unchanged, but the or- 
bital parameters of each binary may be perturbed, possibly 
strongly. Alternatively an exchange encounter may occur, 
where components of the two binaries are exchanged, pro- 
ducing two new binaries. The trajectories of the four stars 
during such an encounter are shown in Figure 1. Another 
possibility is that one binary may be broken up resulting in 
the ejection of two single stars (as illustrated in Figure 2), 
or one star may be ejected leaving the three remaining stars 
in a triple. The formation of a triple-star system is shown in 
Figure 3. Since we only consider hard encounters where the 
total centre-of-mass energy is negative, it is not possible for 
both binaries to be broken up by an encounter. 

Table 2 shows the number of different outcomes for each 
of the runs discussed here. As expected, the effectively larger 
maximum pericentres for encounters with a\/a,2 ^> 1 leads 
to proportionally more flybys in those runs. 

The outcome of each encounter was analysed after the 
run was complete to determine the final state of the bina- 
ries. We do the final analysis of outcome after the runs are 
terminated rather than "on the fly" . It is necessary to do the 
analysis after the runs, because the termination conditions 
are not exact and may lead to excessively long integrations 
for individual encounters, to carry out the simulations in a 
reasonably short time some fraction of unresolved encoun- 
ters must be accepted. The pair of particles with the highest 
specific binding energy was determined; this pair is assumed 
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Table 2. The frequency of the various outcomes. 



Figure 2. An example of a breakup of a binary during a resonant 
encounter between two binaries of unequal mass and semi-major 
axis. 

Stable Triple 



T ' ' ' 1 ' ' 1- 

M 13 =1.0, M 24 =0.5 
a,=3.0, a 2 =0.5, v,/v =l/8-l/4 




Figure 3. A figure showing the ejection of a single star and the 
formation of a triple-star system. As in previous figures the masses 
and initial semi-major axis are unequal. 
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to be bound and its orbital elements were solved for giving a 
new semi-major axis, 03 and eccentricity, e 3 . As the total en- 
ergy was negative, at least one bound pair always exists. The 
other two particles were then considered. If their separation 
from the most tightly bound pair was less than 2az the stars 
were still considered to be strongly interacting and thus the 
encounter was unresolved. If the pair separation was larger, 
the orbital elements of the second pair were solved for, ne- 
glecting any interaction with the first pair. If the resulting 
eccentricity, e4, was less than 1, we concluded two bound 
binaries existed and we solved for the second semi-major 
axis, a a- 

If the separation of the two new binaries was increasing, 
and the kinetic energy of the centres-of-masses of the bina- 
ries exceeded the binary-binary binding energy, we assumed 
the binaries would recede to infinity and we had a flyby or 
exchange, depending on the membership of the respective 
binaries. Else we considered the encounter unresolved. Note 
that this is an approximate condition, unlike the analogous 
case for three body scattering; in practise the separations 
are sufficiently large to ensure that few if any of these cases 
were misclassified as resolved. 

If the second pair of particles was not bound to each 
other, we considered the particle furthest from the most 
tightly bound pair. If this particle was receding from the re- 
maining binary, and its kinetic energy exceeded its binding 
energy relative to the other three particles, it was considered 
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to be escaping; if not, the encounter was still in resonant in- 
terplay and was unresolved. Given that the furthest particle 
was escaping, we now considered the remaining particle. We 
solved for the orbital elements of the third particle about 
the new binary, treating the binary as a point mass at its 
centre-of-mass. If the resulting eccentricity, e t > 1 the third 
star was unbound and we concluded that we had a binary 
and two single stars. If et < 1 the third star was bound 
to the binary and we solved for its semi-major axis rela- 
tive to the binary centre-of-mass, at- To determine whether 
the resulting triple was stable, we used the simple criterion, 
3 x 03(1 + 63) < at(l — et) (Harrington 1975, but see Kiselva 
et aJ. 1994). If this criterion was satisfied, we assumed the 
triple was stable for timescales much larger than the orbital 
timescales, while if it failed we assumed the triple was unsta- 
ble. Unstable triples are expected to decay to a binary and 
a single star, with the new binary semi-major axis drawn 
from a distribution approximately equal to that for reso- 
nant binary-single scattering at the same total energy and 
angular momentum. The approximate nature of the triple 
stability criterion is not a concern at this level of approx- 
imation; relatively few cases are close to being marginally 
stable and if the onset of instability is slow then for triples 
in globular clusters, encounters from field stars are likely to 
determine the future dynamical evolution of the system. A 
similar argument holds for binaries where the energy trans- 
fer was to the binaries and the analysis suggests the stars 
will not escape to infinity. For bound binary trajectories with 
apocentres 3> 03,4, orbital periods become long and in real 
clusters field star perturbations on the quadruple start to 
become significant, the system is no longer dynamically iso- 
lated and for us to continue the integration further becomes 
irrelevant physically. 



3 DYNAMICS AND CROSS-SECTIONS 

For comparison with binary-single calculations and in order 
to calculate physical interaction rates, it is useful to calculate 
(normalised) cross-sections for a process to take place during 
an encounter. This is particularly important for comparing 
runs done with different ai,2, Mi and with different b max . 

We define the cross-section for a process, X, to be 
ax = f{X)i:b'i nax , where f(X) is the fraction of the total 
encounters in that run where X occurred (Hut & Bahcall 
1983, Sigurdsson & Phinney 1993). We define a normalised 
cross-section, bx by 



ox 



ax 



■x{a\ + a\ 



(4) 



This can be compared with the similar definition for binary- 
single scattering (Sigurdsson & Phinney 1993). The normal- 
isation to the binaries geometric cross-section reduces to the 
binary-single case as a 2 — > while the velocity ratio corrects 
for gravitational focusing. However, because of the different 
scaling with a,i for p, the resulting rates are not as simple as 
for the binary-single case. 

A quantity of particular interest is the rate at which a 
given process X occurs in a cluster (core) of density npc~ 3 
and velocity dispersion v s . If there are binaries (of the ap- 
propriate semi-major axis) in the cluster (core), constitut- 
ing some fraction f b of the total number of stars, then we 



expect binaries to encounter each other at some character- 
istic rate, Rx, with corresponding mean time between X 
occurring, Tx = Rx 1 - For a cluster with 1-D line-of-sight 
dispersion v' s , we expect the mean binary encounter veloc- 
ity Vac ~ v s ~ V3v' s , as the binaries encounter each other 
with the full 3-D relative velocity drawn from the under- 
lying velocity distribution, correcting for equipartition. De- 
projection of the low dispersion foreground and background 
contamination, and the weighing of gravitational focusing 
provides additional corrections of order 10% (see Sigurds- 
son & Phinney [1995] for discussion). 

The rate for a process to take place for a given binary 
moving in the cluster (core) is simply 



Rx = (fb-naxv). 



(5) 



It is important in using this approximation to distinguish 
the global binary fraction, that is the fraction of all stars 
that are binaries, from the local binary fraction with the 
relevant range of masses and semi-major axis. Relaxation in 
a cluster may lead to a the fraction of binaries with massive 
primaries in the core being higher than the global binary 
fraction; while dynamical evolution will decrease the core 
binary fraction (see Sigurdsson & Phinney 1995, Hut et al. 
1992). A useful first approximation is to assume a global 
binary fraction of 0.5, with a uniform distribution in semi- 
major axis. Then per decade in semi-major axis, 10% of 
the stars are binaries with a semi-major axis in that range. 
Solving for Tx substituting ox and eliminating v\ we find 



T x = 



(I/Woo)" 1 



, nf b n(ai + a 2 2 )g{M) 



M 1 M 2 _|_ M3M4 



Ox 



(6) 



where g(M) = GM T /(Mi + M 2 )(M 3 + M 4 ). It is useful to 
scale Tx to the cluster parameters of interest. Writing n = 
10 4 ri4 pc~ 3 , = 10«io kms -1 and defining a 1>2 = a lt2 /AU 
and Mi = Mi /M , we can write 



T x = 



1.5 x 10 1 " 

f b g(M)n 4 



((al+al)(E 12 +E 34 )) 7 — \ '^,(7) 

V / \ Voo I Ox 



where we have defined E\ 2 — M\M 2 ja\ and E;n — 
MzM^ja 2 with the masses in solar masses and the semi- 
major axis in AU as before. 

Table 3 shows the normalised cross-section for ex- 
changes, breakups and formation of stable triples for the 
different runs. 

An interesting result is the relatively low cross-section 
for exchange, with a significantly higher cross-section for 
formation of stable triples at all semi-major axis ratios con- 
sidered. As expected the relative cross-section for breakup 
of one of the binaries is large and increases both with Voo/vc 
and 02/01. The fraction of unresolved encounters is small, 
except for the hardest set of runs with ai — a 2 . We don't 
expect the unresolved encounters to contribute dispropor- 
tionately to the close approaches; they mostly consist of 
resonances with one or two stars well separated from the 
other stars, but with insufficient energy to reach infinity; the 
other unresolved encounters consist of unstable triples. We 
expect any "memory" of the initial conditions to have been 
forgotten during the resonances, and the resonances are re- 
solved with a distribution exactly similar to the resonances 
already resolved in that run, providing a few % correction 
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I i ' ' 

M, = 1.0, M 2 = 1.0, M 3 = 1.0, M 4 = 1.0 
number of runs = 4000 
max no. of steps = 500000 
. 0.125 < v/v < 0.25 
P„/a = 5.0/v + 0.6(1.0 + e)v/v c 
e lz = 0.0 

a,' = 1.0, a. = 1.0 . 



" - r 



fit 



Table 3. Cross-sections for outcomes. 



Run a (outcome) 





Exchange 


Breakup 


Triple 


lOd 


0.12 


1.2 


0.40 


lid 


0.23 


2.0 


0.75 


12d 


0.59 


3.5 


1.1 


13d 


0.069 


0.58 


0.24 


20d 


0.14 


1.2 


0.51 


21d 


0.30 


2.3 


0.80 


22d 


0.56 


4.5 


1.0 


lOr 


0.014 


0.14 


0.17 


llr 


0.029 


0.41 


0.41 


12r 


0.044 


1.10 


0.38 


2 Or 


0.012 


0.23 


0.27 


21r 


0.027 


0.42 


0.49 


22r 


0.046 


1.25 


0.24 


30r 


0.049 


0.66 


0.30 


31r 


0.11 


1.2 


0.46 




n 90 

U.iU 


9 9 


u.oo 


40r 


0.053 


0.75 


0.36 


41r 


0.10 


1.3 


0.61 


42r 


0.21 


2.7 


0.83 


50r 


0.0013 


0.011 


0.049 


51r 


0.00086 


0.067 


0.070 


52r 


0.0026 


0.26 


0.014 


60r 





0.010 


0.060 


61r 


0.00065 


0.075 


0.082 


62r 


0.0026 


0.28 


0.022 



to the cumulative cross-sections for close approach at small 
Tmin and adding to the cross-sections for different outcomes 
in the same proportion as the resolved encounters (Heggie 
1988, Sigurdsson & Phinney 1993). Note that for the high 
cii/a? runs we find a significant fraction of breakups of the 
wider binary despite the larger b max at constant C. With 
ai this large, these distant encounters can still significantly 
perturb the wider binary. 



3.1 Close Approach 

A primary goal of this work is to estimate the rate of tidal in- 
teractions during binary-binary encounters in globular clus- 
ters, specifically encounters by hard binaries in the clus- 
ter cores. In order to do that, it is desirable to calculate 
the cumulative (normalised) cross-section for close approach 
between any pair of stars during an encounter, averaged 
over all of the binary phase space. Examples of these cross- 
sections are shown in Figures 4-8, for different mass ratios, 
semi-major axes and velocities. It is useful to fit the (nor- 
malised) cross-sections with a simple broken power-law, as 
shown in Figure 4 (see Hut & Bahcall 1983, Sigurdsson & 
Phinney 1993 for comparable fits for binary-single scatter- 
ing). 

To produce a power law fit, we assume the cumulative 
cross-section, 



lo- 4 io- 3 0.01 0.1 1 

minimum (r ij )/a 

Figure 4. The cumulative cross— section for pairwise close ap- 
proach, m, i = 1 — A,j = i — 4, for particles i,j. The plot shown 
is for equal mass binaries with equal semi-major axis. The sep- 
aration is scaled to the geometric mean of the semi-major axis, 
a = (aic^) 1 ' 2 , and the cross— section normalised to the geomet- 
ric cross-section of the binaries and corrected for gravitational 
focusing. The solid lines show the piecewise power law fit to the 
cross— section, with the dotted line showing the formal breaking 
point between the two fits, tv The error bars are the standard 
error due to finite sampling per separation bin. 

5-(min{r i3 } < r min ) = ■ ( 8 ) 

In the case of binary-single scattering there is only one 
length scale, a, to normalise the "closeness" of approach. 
Here we have two length scales, ai,02, and the normali- 
sation is more ambiguous. We find it useful to scale to the 
geometric mean a = v / aia2 and this is what is plotted in Fig- 
ures 4-8 and is the normalisation used for calculating rates 
later. The choice of a is somewhat arbitrary for oi 7^ <Z2 and 
comparisons of cross-section calculated by different authors 
must be done with care. 

To fit for o"o,7, we applied a piecewise (unweighted) 
least squares fit to the log of the cumulative cross-section 
vs. log Tij. The resulting fits are shown in Table 4. We de- 
fine an additional parameter 77,, the separation at which we 
switch from one power law to the other. For min m < r&, 
(Jo = ffi and 7 = 7!, whereas for min m > rt, 02,72 should 
be used respectively, rt is plotted as a dotted line in Fig- 
ure 4 for illustrative purposes. Formally ri is defined to be 
where the fits cross; in practice this is close to the "knee" 
of the cumulative cross-section curve as desired. Because 
we oversample to large impact parameters, we discard the 
outermost point (nj = 02) in the cumulative cross-section 
when producing the fit. It is necessary to use large enough 
impact parameters to be certain the cross-section for close 
approaches has converged, but for our purposes we are not 
interested in the weak perturbations to the semi-major axis 
caused by the widest flybys. For the high velocity runs with 
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: M, = 1.0, M 2 = 1.0, M 3 = 1.0, M 4 = 1.0 

number of runs = 4000 
- max no. of steps = 500000 
_ 0.5 < v/v c < 1.0 

p„/a = 5.0/v c + 0.6(1.0 + e)v/v c 

e la = 0.0 

a,' = 1.0, a. = 1.0 



lo- 4 io- 3 0.01 

minimum (r^/a 



0.1 



M, = 1.0, M 2 = 0.5, M 3 = 1.0, M 4 = 
number of runs = 4000 
max no. of steps = 500000 
- 0.125 < v/v c < 0.25 

P„/a = 5.0/v + 0.6(1.0 + e)v/v c 

e lz = 0.0 

a,' = 2.0, a. = 0.5 



0.5 



ii 1 



T- 



io- 4 io-= 0.01 

minimum (r^/a 



0.1 



Figure 5. The cumulative cross— section for close approach for the 
high velocity run for equal mass and semi-major axis binaries. 



-TTTT, r " 

: M, = 1.0, M 2 = 0.5, M 3 = 1.0, M 4 = 0.5 

number of runs = 4000 

max no. of steps = 500000 
_ 0.125 < v/v c < 0.25 

p„/a = 5.0/v c + 0.6(1.0 + e)v/v c 

e lz = 0.0 

a,' = 1.0, a_ = 1.0 



III 



io- 4 io- 3 0.01 

minimum (r^/a 



0.1 



Figure 6. The cumulative cross— section for close approach for 
unequal mass and equal semi-major axis binaries. 



Figure 7. The cumulative cross-section for close approach for 
unequal mass and unequal semi-major axis binaries. There is a 
sharp turnup in cross-section for rij = a,2, this occurs because 
of the large number of encounters producing weak perturbations 
to the tighter binary. The discrete binning in rij does not resolve 
the turnup in cross-section. 





M, = 1.0, M 2 = 1.0, M 3 = 1.0, M 4 = 
number of runs = 4000 
max no. of steps — 500000 
. 0.125 < v/v < 0.25 
P„/a = 5.0/v c + 0.6(1.0 + e)v/v c 
e lz = 0.0 

a,' = 4.0, a„ = 0.25 



1.0 



.I'" 



io- 



10- 3 0.01 
minimum (r Vl )/a 



0.1 



Figure 8. The cumulative cross-section for close approach for 
equal mass and unequal semi-major axis binaries. 
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Table 4. The power law fits for cr(r m i n ). 



Run 


5"i 


02 


7i 


72 






n> 


lOd 


130 


5.0 


1.0 


0.24 


1.6 


X 


10~ 2 


lid 


160 


4.8 


1.1 


0.27 


1.8 


X 


10~ 2 


12d 


210 


5.7 


1.3 


0.41 


1.7 


X 


10~ 2 


13d 


160 


5.2 


1.1 


0.24 


1.6 


X 


10~ 2 


20d 


150 


6.0 


1.1 


0.32 


1.8 


X 


10~ 2 


21d 


230 


6.4 


1.3 


0.39 


2.0 


X 


10~ 2 


22d 


120 


6.8 


1.2 


0.46 


2.0 


X 


10~ 2 


lOr 


28 


0.61 


1.2 


0.42 


7.4 


X 


10~ 3 


llr 


9.5 


0.73 


0.98 


0.39 


1.3 


X 


10~ 2 


12r 


21 


0.49 


1.3 


0.41 


1.5 


X 


10~ 2 


20r 


5.8 


0.9 


0.86 


0.35 


2.6 


X 


10~ 2 


21r 


4.6 


0.86 


0.86 


0.46 


1.5 


X 


10~ 2 


22r 


4.0 


0.61 


0.98 


0.49 


2.2 


X 


10~ 2 


30r 


80 


2.2 


1.1 


0.26 


1.4 


X 


10~ 2 


31r 


40 


2.1 


1.0 


0.31 


1.4 


X 


10~ 2 


32r 


20 


2.2 


1.0 


0.47 


1.6 


X 


10~ 2 


40r 


49 


2.7 


1.1 


0.35 


2.1 


X 


10~ 2 


41r 


34 


2.6 


1.0 


0.37 


1.7 


X 


10~ 2 


42r 


96 


2.8 


1.3 


0.49 


1.3 


X 


10~ 2 


50r 


0.27 


0.031 


0.9 


0.43 


1.0 


X 


10~ 2 


51r 


15 


0.019 


1.8 


0.42 


8.0 


X 


10~ 3 


52r 


0.038 




0.95 










60r 


0.41 


0.027 


1.1 


0.43 


1.7 


X 


10~ 2 


61r 


0.11 


0.027 


0.94 


0.55 


2.7 


X 


10~ 2 


62r 


0.048 




0.89 











ai = 16a2 a single power law fit was adequate for all m as 
indicated in the table. 

The error-bars shown in Figures 4-8 show the standard 
deviation due to sampling noise from the finite number of 
encounters per bin. As the cross-section is cumulative the 
errors are systematically smaller at larger m , hence it is not 
appropriate to weight the least squares fit by the errors. 

As all the runs presented here are for ei,2 = there 
might be some concern that for binaries with a thermal 
eccentricity distribution, as might be expected in globular 
clusters, the cross-section for close approach could be sub- 
stantially larger than inferred from these simulations. To 
explore this we performed some runs with ei,2 drawn from a 
thermal probability distribution, P(e;) = 2e», and compared 
the ratio of cumulative cross-section for close approach for 
the eccentric and zero eccentricity binaries. 

The result is shown in Figure 9. As can be seen from 
the figure, the ratio is approximately 1.7 = 1+ < e >, over 
a large range in Vijja. The ratio rises at small rij /a, and re- 
sults from other runs with ai / ai and ei,2 thermal suggest 
this is significant. The ratio approaches unity at rij /a = 1 by 
definition. We conclude that allowing for an eccentricity dis- 
tribution increases the cross-section for collisions by a factor 
of ~ 1.7, except for Vij/a ~ 10~ 4 for which the correction 
is somewhat larger. The same effect is seen for 02 7^ ai but 
somewhat more pronounced. A possible concern is that this 
not a real dynamical effect, but simply the high eccentricity 
end of 122(1 — e.2) periastron passage. To check that this was 



o 

II 



M 13 = 1, M 24 = 1 
a, = 1, a, = 1 
- 0.135 < v/v < 0.25 



— 1 



10- 4 10- 3 0.01 0.1 1 

minimum (rv/a) 

Figure 9. The ratio of cumulative cross— section for close ap- 
proach for the for equal mass, equal semi-major axis binaries for 
zero eccentricity and for P(e) = 2e runs. The errorbars show the 
standard error due to finite sampling. 



not the case we made two checks: we truncated the eccen- 
tricity distribution at ei,2 = 0.98 and checked the increase 
in relative cross-section was still present, as it was, and we 
compared the initial periastron separation for each binary 
with min rij/a. The closest approaches were predominantly 
due to binaries with ei,2 ~ 0.5, not encounters with very 
high initial binary eccentricity. 



3.2 Hardening of Binaries by Fly-by 
Encounters 

We now consider the effects on the binaries of fly-by encoun- 
ters. Such events are common compared to other processes 
such as exchanges, or the formation of triples as shown in 
Table 2. In Figures 10 and 12, we plot the distribution of 
semi- major axes a[, and a' 2 of the two binaries after a fly- 
by encounter for 01 = 02 and 01 = 4a2 respectively. As 
expected the distribution is symmetric for the equal semi- 
major axis binaries. The empty region in the upper right 
hand corner is forced by energy conservation, while angular 
momentum conservation prevents both binaries from hard- 
ening a lot simultaneously. For ai = 4a2 most of the encoun- 
ters produce a very weak perturbation in ai compared with 
01, as expected. 

Figures 11 and 13 show the differential cross-sections 
for change in semi-major axis (cf. Davies et al. 1994, Sig- 
urdsson & Phinney 1993, Hut & Bahcall 1983). For binary- 
single scattering, we expect da/dA oc (1 + A) -4 ' 5 for strong 
encounters (Heggie 1975), where A = SEb/Et, is the frac- 
tional change in binary binding energy. In semi-major axis 
space we expect da/dai oc af 5 for strong encounters. The 
dotted line in Figures 11 and 13 shows a af 5 power law, 
the differential cross-section shows an approximate fit to 
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10 3 



10 a 



b 

T3 10 



0.1 



0.1 



a, = 1, a 3 = 1 
0.125 < v,/v < 0.35 



10 

log(a,/a) 



10 2 



10 3 



Figure 10. Plot of log CI3 against log 04, for flybys for equal mass 
and initial semi-major axis binaries. There is a "forbidden" re- 
gion excluded by energy conservation, bounded by the dotted line 
in the figure. The distribution in 03,04 is roughly symmetric as 
expected. The changes in semi-major axis included 775 encoun- 
ters where 03 < a± and < 04, 139 encounters where both 
binaries increased semi— major axis, and 549 and 591 encounters 
where one binary increased its semi— major axis and the other 
decreased, respectively. 



Figure 11. Distribution of final semi— major axis. The solid his- 
togram show the distribution for 03, the dashed line shows the 
distribution for 04. The cross-section is unnormalised. The dashed 
line shows the distribution in ai expected from Heggie's law for 
binary— single scattering, da/dai ex. a?' 5 . The distribution shown 
in the figure is for equal mass binaries and cii = 0,2- As expected, 
the distributions in 03,04 arc equal to within poisson noise. The 
spike in the distribution at 034 = 1 is due to the large number 
of wide encounters producing small changes in semi-major axis. 



Heggie's law over an interesting range in 8a/ a for a suf- 
ficiently large 5a. For weak perturbations, the differential 
cross-section diverges with b max whereas for very strong en- 
counters the falloff in cross-section is somewhat steeper than 
expected from Heggie's law. It is interesting to note that res- 
onant flybys do not lead to large changes in binding energy 
compared to binary-single encounters. For other classes of 
outcome we defer discussion of the changes in energy to a 
later paper. 



3.3 Application of Cross Sections to Compute 
Rates in Globular Clusters 

We now apply the cross-sections computed above to models 
of globular clusters. The calculated rates are listed with the 
cluster parameters in Table 5. Substituting equation 8 into 
equation 7, and approximating (l/uoo) -1 = v s we get the 
rate of collisions 



1.5 x 10 1U 

f b g(M)mv : 



(ial + ai){E 12 + 



£34)) 



_1 (0102 



,7/2 



cro(r„ 



■y,(9) 



where 7 = 71 , ao = o\ for r min < 77, , and 7 = 72 , <7o = ai 
for Tram > TV Table 5 shows a representative set of Tc for 
a range of cluster parameters, and for binaries containing 
two O.7M main-sequence stars and for binaries containing 
a 1.4Mq compact star and a 0.7Mq main-sequence star re- 
spectively. The binary separation assumed is shown for 01 in 
AU, with ai/a2 then as given in column 5. For the calcula- 
tion we choose r m i n = ftR*, where R, is the stellar radius, 




Figure 12. Plot of 03 against 04, for flybys for equal mass and un- 
equal initial semi-major axis binaries. The dotted line shows the 
"forbidden" region excluded by energy conservation. The changes 
in semi— major axis included 973 encounters where 03 < ai and 
02 < 04, 770 encounters where both binaries increased semi- 
major axis, 648 encounters where 03 > a\ and 04 < 02, and 1094 
encounters where 03 < ai and 04 > 03. 
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Table 5. The timcscale for tidal interaction or stellar collision, 
Tc (in years), for a binary with semi-major axis 01 encountering a 
binary with semi-major axis 02 in various globular cluster models, 
assuming f b = 1 in the core for the semi-major axis used. a\ is 
shown in AU, and M* is the mass of the more massive star in 
each binary in solar masses. v s is given in kms -1 . 
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Figure 13. Distribution of final semi— major axis. The solid his- 
togram show the distribution for 0,3, the dashed line shows the 
distribution for 0,4. The cross-section is unormalised. The dashed 
line shows the distribution in expected from Heggie's law for 
binary— single scattering, da/dai ex. af' 5 - The distribution shown 
here is for equal mass binaries and a\ = 4a 2 , as expected the 
distributions in 03,04 arc different, 04 showing a much narrower 
spread from its initial value. Rather surprisingly, Heggie's law 
still fits well for strong perturbations to the tighter binary, for 
04 ~ l/2a2. The spike in the distribution at 03 =2,04 = 0.5 
is due to the large number of wide encounters producing small 
changes in semi-major axis. 



and ft ~ 3 for main-sequence stars, with ftR* being the 
separation at which tidal effects become significant in these 
encounters (Press & Teukolsky 1977). As before, the masses 
are in solar masses and radii in AU. The timescales in Table 
5 are calculated assuming /(, = 1, that is all the stars in the 
cluster (core) are binaries with the appropriate mass and 
semi-major axis. In practise fb < 1, possibly much less, and 
the total collision rate represents an average over all binary 
masses and semi-major axes. 

A typical globular cluster core will contain about 



10 stars. Here we follow observational conven- 



tions and count a binary as a single "star" when figuring 
N*. Approximating the total collision rate as coming from 
a constant density core and neglecting contributions due to 
binary interactions outside the core, an approximation for 
the total number of collision products observed in a cluster 
is given by 



Nc = fefdf b 2 N^ y 



(10) 



where r is the characteristic lifetime over which the collision 
product (eg. a blue straggler) is observable. Here we sepa- 
rate the fractional binary density in the core, f d and the 
binary fraction in the appropriate semi-major axis range, 
fb- We also allow for a factor f e to correct for the eccen- 
tricity distribution. For most cases / e ~ 1.7, except noting 
that for R„ £ 10~ 3 a 2 and a 2 <§C 01, f e ~ 10. The true 
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rate is the integrated rate over T c (Mi,cn) given f b (Mi,ai); 
however the uncertainty in the binary population distribu- 
tion is very large and the systematic uncertainties in the 
integrand make a formal integral meaningless. An approxi- 
mation to the true collision rate can be made by assuming 
that 0(6%) of the stars are near the turnoff, and another 
0(2 — 4%) is in evolved remnants. Mass-segregation then 
increases the fractional density of these more massive stars 
in the core further. As a simple approximation we can also 
assume that the binary fraction, f b , is about 0.1 per decade 
in semi-major axis independent of binary mass, and that 
the total initial binary population spans about 5 decades in 
semi-major axis, a, ~ 10~ 2 — 10 3 AU. 

In an old globular cluster, the core binary population is 
dynamically modified by mass segregation and subsequent 
interactions. Hardening and breakup reduce the number of 
wide binaries but increase the fraction of binaries with semi- 
major axis, ~ a c , such that Tf}(a c ) ~ tv, where r r is the 
cluster relaxation time and Tn(a c ) is the characteristic time 
scale for a binary with semi-major axis a c to change its 
semi-major axis by order a c /2. The situation is further com- 
plicated by dynamical recoil during strong binary-single and 
binary-binary interactions, which tends to remove hardened 
binaries from the cluster core to the cluster halo where in- 
teraction timescales are long (see eg. Sigurdsson & Phinney 
1995). Binaries ejected to the cluster halo then return to the 
cluster core on a relaxation timescale. 

As can be seen from Table 5, the collision timescale 
for a binary of some semi-major axis a\ interacting with 
a binary of semi-major axis ai < ai does not vary much 
from a\ = ai to ai ~ 10a2, but then becomes longer for 
ai 3> 10<Z2- Thus we can use a meaningful average colli- 
sion timescale for binaries with semi-major axis in the same 
decade span, and neglect collisional interaction for binaries 
with semi-major axis in different decade spans as being neg- 
ligible by comparison. Consider three different clusters, with 
core densities, n, of 10 2 ,10 4 and 10 6 pc~ 3 respectively. As- 
suming typical concentration parameters of Wo ~ 6, 9, 12 for 
multi-mass Michie-King models and a Salpeter initial mass- 
function, mass-segregation increases the fractional density 
of 0.7 + O.7M0 binaries of all semi-major axes from 0.06, 
to 0.08,0.15 and 0.26 respectively, neglecting dynamical re- 
coil. Typical r r are 10 10 , 10 9 and 10 8 years respectively, 
and we expect those main-sequence binaries dominating the 
interaction rate to have semi-major axes of 10 — 100AU, 
0.1 — 1AU and < 0.1 AU respectively. Assigning a binary 
fraction of f b (a t = 10 - 100 AU) « 0.2 for the low den- 
sity cluster, f d = 0.08 and f e = 10, we find the expected 
number of currently observable main-sequence collisions to 
be, N BS (n = 10 2 pc~ 3 ) = 4(/ 6 /0.2) 2 per 10 4 core stars for 
tbs — 5 x 10 9 years, and taking T c w 10 11 years as indicated 
in Table 5. For the medium density cluster, we similarly find 
N BS (n = 10 4 pc~ 3 ) = 40 for f b ( ai = 0.1 - 1AU) = 0.2, 
fd = 0.15 and f e = 1.7. A higher f b may be appropriate as 
the interaction time scale is short enough for a substantial 
fraction of binaries to be hardened from larger semi-major 
axes to the optimum range for collisions. For the densest 
cluster, we expect fd to be smaller than mass-segregation 
would indicate as dynamical recoil and breakup are likely 
to have been significant over the cluster history, so we take 
f e = 1.7, f b (oi < 0.1AU) = 0.1 and f d = 0.1, even though 
mass-segregation would suggest f d ~ 0.2 — 0.3. Using these 



values we find JVas(n = 10 pc ) = 500. However, the 
true number may be smaller still as exchanges and collisions 
may have reduced the fraction of core stars near the main- 
sequence turnoff relative to the evolved remnant stars (see 
Sigurdsson & Phinney 1995). 

It is interesting to compare, for binaries in the same 
semi-major axis range, the timescales to breakup due to 
binary-binary interactions with the timescale to collision. 
Using the results in Table 3, and binary parameters as 
above, we find T brea k up (n — 10 2 pc~ 3 ) = 3 x 10 10 years, 
T breakup (n = 10 4 pc~ 3 ) = 4 x 10 1() years and T breakup (n = 
10 6 pc" 3 ) = 3 x 10 9 years. For the lowest density cluster, 
a binary is several times more likely to be broken up by 
a binary-binary interaction, than be involved in a stellar 
collision, whereas for the medium density cluster breakup 
and collision are about equally likely, and the chance of a 
breakup in a binary-binary collision in the densest cluster 
is somewhat less than the probability of a stellar collision. 



3.4 Comparison with Previous Results 

NGC 5053 is a very low density (n = 8 pc -3 ), low dispersion 
globular, with 24 observed candidate bue stragglers (Nemec 
& Cohen 1989). Some controversy exists in the literature 
over whether the blue stragglers in this cluster must be 
due to merger through internal evolution of initially tight 
primordial binaries (Hills & Day 1976) or whether binary- 
binary collisions may have produced a significant fraction 
of the observed blue straggler stars (Leonard & Fahlman 
1991). 

We find that we cannot produce the 24 blue stragglers 
observed in NGC 5053 even with f b = 1 for binaries with 
oi ~Z 100AU; rather we find, T c « 2 x 10 12 years, giving 
Nbs = 2 - 3(/ e /10)(/ d /0.1)(/f,/0.2) 2 for a total core popu- 
lation of 40,000 stars. The encounter rate is dominated by 
binaries with semi-major axes <; 100 AU, and it is unlikely 
fb ^> 0.2 for this semi-major axis range, as many binaries 
must be much tighter than this. It is very unlikely that mass 
segregation has increased the fraction of binaries with pri- 
maries of near turnoff mass to 3> 0.1 in this cluster; the 
relaxation time scale is very long and the concentration low 
(c = 0.75). 

It is not necessary to produce all the blue strag- 
glers through binary-binary collisions, as some undoubtedly 
formed through spiral-in of tight primordial binaries (see eg. 
Livio 1993). It is possible that the fraction of primordial bi- 
naries with initial semi-major axis small enough for them 
to merge in a globular cluster lifetime is a function of the 
globular cluster initial density and dispersion. It is also pos- 
sible that the core density of low density globular clusters 
has decreased on timescales of few x 10 years through tidal 
shocking by the galaxy, and that the collision rates were 
thus higher in the past. Note that for every blue straggler 
produced by collisions during binary-binary encounters in 
these globular clusters we expect <; 10 binaries to have been 
dynamically broken up, and thus the current binary frac- 
tion in the right semi-major axis range would be smaller 
now than in the past. 

In the case of NGC 5053, we might require the core 
density to have been an order of magnitude higher within 
the lifetime of the currently existent blue stragglers, and a 
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large fraction (<; 0.3) of the core turnoff stars to have been 
in wide (a ~ 100AU) binaries. We do expect the binary 
fraction to have decreased with time, but we still expect the 
present core binary fraction to be ~ 0.25 if binary-binary 
collisions are to account for the blue stragglers. This would 
be apparent in high precision photometry of NGC 5053 as 
a prominent second main-sequence (Romani & Weinberg 
1991). Alternatively low density clusters like NGC 5053 may 
form a significant number of wide triples, and collisions are 
frequent during triple encounters. 

It is instructive to compare the collision rate due to 
binary-binary encounters with the collision rate due to 
binary-single encounters. For the sample clusters discussed 
above, we find binaries with those semi-major axes have 
collision timescales with a single main-sequence star of 
4.6 x 10 11 , 5.4 x 10 9 and 1.8 x 10 s years respectively, and 
the expected number of blue stragglers formed is given by 
Nbs = fdfbiX — fb)f*N*r/T c , where /» is the fractional 
density of main-sequence stars at the turnoff in the core, 
corrected for mass segregation (0.13,0.2,0.2 respectively). 
Using the same f b as above, we find the expected number of 
blue stragglers to be — 1,50 and 1100 respectively. That 
is, in the low density clusters binary-single collisions rates 
are comparable to binary-binary collisions, and the binary- 
singles dominate in the denser clusters. As not that many 
blue stragglers are observed in globular clusters one might 
infer the binary fraction per decade in semi-major axis is less 
than the 0.1 used here and the global binary fraction some- 
what less than 0.5, which is consistent with observational 
estimates (Pryor et al. 1989, Yan & Mateo 1994) while still 
allowing a sufficiently high collision rate to produce the blue 
stragglers and other stellar exotica observed. 



4 CONCLUSIONS 

Some care must be taken in considering the effectiveness of 
binary-binary collisions in globular clusters. The global bi- 
nary fraction at zero age in clusters is probably 0.5 — 1.0, 
comparable with that seen in the field. However, this in- 
cludes binaries from near contact, a; ~ 0.01AU, to ex- 
tremely wide binaries, a; 3> W^AU. The former do not 
interact on short enough timescales to be of interest, ex- 
cept in core-collapsed clusters, and will in due course merge 
through their internal evolution; the latter are soft and have 
high encounter rates in all except the very lowest density 
clusters, and are broken up in a few dynamical timescales. 
The global binary fraction as a function of <Zi seems to have 
an initial distribution of approximately 0.1 per decade in a;, 
and that is the approximation we use above. However, the 
core population of binaries is modified by several processes, 
including mass-segregation, dynamical recoil, exchange and 
breakup. As a result the fraction of binaries in the core, per 
decade in Oj and at different masses varies with time and 
cluster parameters. In calculating the expected number of 
blue stragglers above we made some effort to correct for the 
dominant processes in the different clusters considered. 

We find collision timescales for plausible binary popu- 
lations comparable to the lifetime of the clusters, and an 
expected number of blue stragglers sufficient to account for 
a large fraction of the low density blue straggler population, 
but overestimating the population in the denser clusters. 



This can be understood in terms of the dynamical evolution 
of the globular cluster binary population, as breakup and 
ejection decreases the core population of binaries. 

We have refrained here from discussing in detail the 
properties of the final state of the binaries. In particular, 
parameters of interest include the final distribution of semi- 
major axis, not just for the fiybys and exchanges, but also 
the breakups and triples, and the resultant cross-sections 
for energy transfer and recoil velocity distribution. Also of 
interest are the eccentricity distributions of the various final 
binary states. Of particular interest to us are the properties 
of the system after it undergoes an inelastic collision. Sim- 
ulations of such collisions have been performed using SPH 
(Davies et al. 1993, 1994, Goodman & Hernquist 1991, Sig- 
urdsson & Hernquist 1993). Approximating the collision as a 
totally inelastic "sticky particle" merger, conserving momen- 
tum but not energy, allows a quick and reasonably accurate 
way of determining the properties of the merged systems, in 
particular whether they form a single merged star, or if the 
merged star is in a binary or even a triple, and if so what the 
orbital parameters and center of mass recoil velocity of the 
system containing the merged star is. An analysis of these 
properties is deferred to a second paper (in preparation). 

It is clear that binary-binary interactions are significant 
for producing stellar exotica through collisions in globular 
cluster cores. Compared to binary-single interactions, the 
rates inferred suggest a modest global binary fraction in the 
cores of the denser clusters, in accord with previous esti- 
mates, with /b ( all <ii) ~ 0.2 and fb(ai) ~ 0.05 per decade in 
di, while in the low density clusters the blue straggler popu- 
lation is consistent with a somewhat higher binary popula- 
tion, with perhaps > 10% of the turnoff mass main-sequence 
stars in the core being in binaries with a; ~ 100AU. Binary- 
binary collisions most likely dominate binary-single colli- 
sions in many low density clusters as suggested by Leonard 
(1989) and may account for a significant fraction of the blue 
stragglers observed. 
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